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Final  Report  W911NF-09-1-0147:  High  performance  nuclear  magnetic 
resonance  imaging  using  magnetic  resonance  force  microscopy 


P.  Chris  Hammel 

Department  of  Physics ,  The  Ohio  State  University 

We  summarize  the  key  advances  resulting  from  our  research  program  supported  by  ARO  award 
W911NF-09-1-0147  between  July  15,  2009  and  July  14,  2013.  We  have  used  ultra-sensitive  Magnetic 
Resonance  Force  Microscopy  detection  Electron  Spin  Resonance,  using  (ESR-MRFM),  to  perform 
spatially  resolved  scanned  probe  studies  of  spin  dynamics  in  nanoscale  ensembles  of  few  electron 
spins  of  varying  size.  Our  research  culminated  in  observation  that  the  distinction  between  two 
channels  of  energy  flow  from  and  within  a  spin  system,  typically  characterized  the  spin-lattice  re¬ 
laxation  time  Ti  and  the  spin-spin  relaxation  time  T2  that  are  rigorously  distinct  in  macroscopic 
systems  become  indistinguishable  in  nanoscale  systems.  This  is  because  spin  transport  mediated 
by  mutual  spin  flip-flop  diffusion — arising  from  a  T2- like  process — in  a  nanoscale  quasi- ID  spin 
nanowire  can  transfer  energy  from  the  measured  spin  system  to  the  thermal  reservoir  outside  the 
measured  volume — behavior  typical  of  Ti  processes.  The  spin  nanowire  we  created  in  diamond 
by  means  of  patterned  nitrogen  ion  implantation  represents  an  intriguing  prototype  element  of  a 
spintronic  device  in  an  attractive  and  widely  used  wide  band  gap  semiconductor.  ESR-MRFM,  a 
scanned  probe  method  capable  of  spatially  resolved  magnetic  resonance  detection,  offers  a  power¬ 
ful  approach  to  understanding  the  spin  properties  and  dynamics  in  the  device-like  environments 
where  their  dynamics  and  lifetime  will  be  influenced  by  nanoscale  confinement.  Our  theoretical 
model  that  explains  the  observed  spin  dynamics  in  confined  nanoscale  spin  structures  are  sup¬ 
ported  by  extensive  numerical  modeling  that  agree  well  with  our  experimental  data.  Using  our 
understanding  of  the  electron  spin  dynamics  in  confined  nanoscale  structures,  we  developed  an 
improved  measurement  protocol,  pi-OSCAR,  that  we  that  we  used  to  demonstrate  spectroscopic 
measurements  on  nanoscale  spin  ensembles  of  few  electron  spins.  Our  demonstration  of  nanoscale, 
spatially  resolved  measurements  of  spin  dynamics  in  quasi  ID  spin  nanowires  in  diamond,  pre¬ 
sented  in  a  manuscript  entitled  “The  effect  of  spin  transport  on  lifetime  in  nanoscale  systems,” 
is  currently  under  review  at  Nature  Nanotechnology.  This  advance  involved  several  advances: 

•  Creation  of  an  integrated  model  of  spin  dynamics  and  spin  transfer  across  interface  in 
nanoscale  spin  ensembles  in  the  presence  of  large  magnetic  field  gradients. 

•  Microscopic  Monte  Carlo  simulation  of  spin-flip  mediated  spin  transport  in  quasi  ID  systems. 

•  Development  of  a  data  deconvolution  protocol  enabling  spectroscopic  studies  of  nanoscale 
spin  ensembles. 

•  Development  and  detailed  analysis  of  a  new  ESR-MRFM  detection  protocol — pi-OSCAR — 
that  improves  control  over  the  size  of  the  measured  spin  ensemble. 

These  achievements  were  enabled  by  improvements  in  ESR-MRFM  technique  and  detec¬ 
tion.  These  included  improvement  of  the  signal-to-noise  ratio  of  the  measurement  appara¬ 
tus,  preliminary  studies  of  the  material  systems  for  the  current  and  next-step  experiments 
with  the  ESR-MRFM,  and  fabrication  of  the  nitrogen  implanted  spin  nanowires;  these  include: 

•  Detection  of  ~  13  electron  spins  with  high  signal-to-noise  ratio  through  improved  MRFM 
cantilever  visibility,  quality  factor  and  microwave  delivery. 

•  Custom  N-implanted  diamond  samples  including  uniformly  doped  samples  and  quasi  ID  spin 
nanowires  were  fabricated  for  ESR-MRFM  imaging. 

•  PI  and  NV  centers  in  nitrogen  implanted  diamond  samples  were  simultaneously  measured 
demonstrating  the  ability  to  detect  and  characterize  defect. 

•  Conventional  EPR  and  ENDOR  characterization  of  PI  centers  in  diamond  for  nuclear  spin 
detection  through  electron  spin  detection. 

•  Studies  of  interfacial  magnetism  in  LaA103/SrTi03(LA0/ST0). 
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I.  INTRODUCTION 

We  summarize  the  key  advances  resulting  from 
our  research  program  supported  by  ARO  award 
W911NF-09-1-0147  between  July  15,  2009  and 
July  14,  2013.  We  have  used  ultra-sensitive  Mag¬ 
netic  Resonance  Force  Microscopy  detection  Elec¬ 
tron  Spin  Resonance,  using  (ESR-MRFM),  to 
perform  spatially  resolved  scanned  probe  stud¬ 
ies  of  spin  dynamics  in  nanoscale  ensembles  of 
few  electron  spins  of  varying  size.  Our  research 
culminated  in  observation  that  the  distinction 
between  two  channels  of  energy  flow  from  and 
within  a  spin  system,  typically  characterized  the 
spin-lattice  relaxation  time  T\  and  the  spin-spin 
relaxation  time  that  are  rigorously  distinct  in 
macroscopic  systems  become  indistinguishable  in 
nanoscale  systems.  This  is  because  spin  trans¬ 
port  mediated  by  mutual  spin  flip-flop  diffusion — 
arising  from  a  TVlike  process — in  a  nanoscale 
quasi- ID  spin  nanowire  can  transfer  energy  from 
the  measured  spin  system  to  the  thermal  reservoir 
outside  the  measured  volume — behavior  typical 
of  Ti  processes.  The  spin  nanowire  we  created 
in  diamond  by  means  of  patterned  nitrogen  ion 
implantation  represents  an  intriguing  prototype 
element  of  a  spintronic  device  in  an  attractive 
and  widely  used  wide  band  gap  semiconductor. 
ESR-MRFM,  a  scanned  probe  method  capable  of 
spatially  resolved  magnetic  resonance  detection, 
offers  a  powerful  approach  to  understanding  the 
spin  properties  and  dynamics  in  the  device-like 
environments  where  their  dynamics  and  lifetime 
will  be  influenced  by  nanoscale  confinement.  Our 
theoretical  model  that  explains  the  observed  spin 
dynamics  in  confined  nanoscale  spin  structure- 
sare  supported  by  extensive  numerical  modeling 
that  agree  well  with  our  experimental  data.  Using 
our  understanding  of  the  electron  spin  dynam¬ 
ics  in  confined  nanoscale  structures,  we  developed 
an  improved  measurement  protocol,  pi-OSCAR, 
that  we  that  we  used  to  demonstrate  spectro¬ 
scopic  measurements  on  nanoscale  spin  ensembles 
of  few  electron  spins.  The  results  of  our  studies 
will  be  presented  in  a  manuscript,  “The  effect  of 
spin  transport  on  lifetime  in  nanoscale  systems,” 
currently  under  review  at  Nature  Nanotechnol¬ 
ogy- 


These  achievements  were  supported,  particu¬ 
larly  in  the  earlier  stages  of  the  project,  by  im¬ 
provements  of  the  ESR-MRFM  technique.  These 
included  improvement  of  the  signal-to-noise  ra¬ 
tio  of  the  measurement  apparatus,  preliminary 
studies  of  the  material  systems  for  the  current 
and  next-step  experiments  with  the  ESR-MRFM, 
and  fabrication  of  the  nitrogen  implanted  spin 
nanowires. 

These  studies  have  advanced  ability  to  study 
spin  properties  of  nanoscale  structures  and  have 
explored  new  approaches  to  spintronics  in  dia¬ 
mond,  a  wide  bandgap  semiconductors  that  could 
be  attractive  for  scientific  and  technological  prob¬ 
lems  of  importance  to  the  DoD,  industry,  and 
society.  This  research  is  based  on  our  unique 
capability  for  ultrasensitive  mechanical  detection 
of  magnetic  resonance.  This  research  has  bene¬ 
fited  from  a  close  interaction  with  our  collabora¬ 
tors  including  Dr.  D.  Rugar  at  IBM  (Almaden, 
CA),  Prof.  R.  Budakian  at  UIUC  (Urbana- 
Champaign,  IL)  and  Prof.  C.  H.  Ahn  at  Yale 
University  (New  Haven,  CT). 

The  most  interesting  result  of  our  studies  is 
the  ESR-MRFM  study  of  spin  dynamics  quasi  ID 
spin  nanowires  in  diamond.  The  measurements 
show  that  the  dynamics  of  the  electron  spins  in  a 
nanoscale  ensemble  are  dominated  by  transfer  of 
the  spin  moment  across  the  spatial  boundaries  of 
the  measurement  volume  via  spin  diffusion  facili¬ 
tated  by  mutual  spin  flip  transitions.  Also  inter¬ 
esting  is  the  fact  that  these  measurements  were 
conducted  by  coupling  spin  magnetic  moments 
to  the  high  gradient  magnetic  field  of  the  ESR- 
MRFM  probe  which  will  suppress  the  mutual  spin 
flip  diffusion.  Our  investigation  of  this  seeming 
contradiction  led  us  to  consider  several  expla¬ 
nations  for  the  observed  effect:  these  included 
the  possibility  that  the  strongly  anisotropic  hy- 
perfine  interaction  of  an  electron  on  PI  center 
with  the  associated  14 N  atom  could  be  respon¬ 
sible  for  overcoming  this  suppression.  To  resolve 
this  question,  we  conducted  an  extensive  numer¬ 
ical  microscopic  Monte  Carlo  simulation  of  the 
spin  diffusion  assisted  spin  transport  in  quasi  ID 
systems  (see  Sec.  Ill)  and  developed  a  method  for 
ESR-MRFM  data  deconvolution  enabling  spec¬ 
troscopic  measurements  on  nanoscale  spin  ensem- 
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bles  (see  Sec.  IV).  Our  investigations  show  that 
the  observed  absence  of  spin  diffusion  suppression 
stems  from  the  inhomogeneous  ESR  line  broad¬ 
ening  of  an  electron  on  PI  center  due  to  the  dis¬ 
tortions  introduced  by  ion  implantation.  These 
findings  are  summarized  in  Sec.  II.  We  have 
also  conducted  a  detailed  analysis  of  a  new  ESR- 
MRFM  spin  manipulation  protocol  that  allows 
improved  control  over  the  measurement  volume 
in  ESR-MRFM  experiment  (see  Sec.  V). 


II.  SPIN  DYNAMICS  IN  NANOSCALE 
SPIN  ENSEMBLES 

Spin  centers  in  diamond — such  as  nitrogen  va¬ 
cancy  (NV)  centers,  nitrogen  (PI)  centers,  or  the 
newly  discovered  ST1  center  [1] — are  of  great  in¬ 
terest  for  spintronics  and  quantum  computation 
and  have  proven  to  be  valuable  platforms  for 
studying  the  role  of  spin  interactions  [2].  Spins 
which  are  well-isolated  from  their  environment 
exhibit  long  coherence  times  [3-5],  yet  spin  in¬ 
teractions  are  necessary  for  entanglement,  trans¬ 
port,  and  read-out  schemes  [6].  Here,  we  directly 
observe  spin  transport  in  dipole-coupled  PI  spins 
in  thermal  equilibrium.  Since  diamond  is  an  in¬ 
sulator,  charge  motion  is  completely  absent  from 
this  spin  transport.  The  transport  is  confined 
to  a  nano-scale  quasi- ID  ‘spin  wire’  where  the 
spin  density  leads  to  appreciable  dipolar  coupling 
and  spin  flip-flop  interactions.  In  samples  with 
long  T\  —  such  as  our  own  —  this  type  of  spin 
transport  could  result  in  highly  efficient,  even  re¬ 
versible  [7],  pure  spin  currents.  Furthermore,  by 
implementing  a  magnetic  resonance  protocol  that 
improves  spatial  resolution,  we  demonstrate  se¬ 
lective  detection  of  different  nuclear  spin  states 
and  measure  the  hyperfine  spectrum  within  the 
spin  wire.  Our  measurements  probe  intrinsic  spin 
dynamics  at  the  nanometer  scale,  providing  crit¬ 
ical  information  for  practical  devices  which  seek 
to  control  spin. 

Spintronics  research  has  largely  focused  on  two 
regimes  —  spin  transport  by  conduction  elec¬ 
trons,  modeled  as  a  continuum  spin  density  [8], 
and  the  isolated  single  spin  “qubit”  [9,  10].  Ap¬ 
propriate  coupling  between  spins  could  allow  for 


efficient  spin  transport,  while  preserving  long  spin 
lifetimes.  For  example,  a  spin  flip-flop  —  the 
simultaneous  flipping  of  a  pair  of  anti-aligned 
neighboring  spins  due  to  their  dipolar  interaction 
—  conserves  magnetization.  Hence,  successive 
flip-flops  can  result  in  a  pure,  though  diffusive, 
spin  current  (i.e.  not  accompanied  by  charge  mo¬ 
tion).  In  our  work,  we  show  evidence  of  pure 
spin  transport  —  an  often  sought-after  goal  for 
spin  transport  devices  —  across  defect  centers  in 
diamond.  Our  measurements  are  performed  on 
an  ensemble  containing  fewer  than  100  net  spins 
under  conditions  where  the  thermally-averaged 
polarization  and  polarization  gradient  are  vanish¬ 
ingly  small.  In  contrast  to  non-equilibrium  exper¬ 
iments,  monitoring  the  ‘spin  noise’  —  stochastic, 
temporary  deviations  of  the  polarization  which 
are  “averaged  out”  in  measurements  of  large  en¬ 
sembles  —  allows  us  to  characterize  the  intrinsic 
spin  dynamics  of  the  ensemble  [ii]. 

For  small  spin  ensembles,  fluctuations  become 
the  key  to  not  only  measuring  spin  diffusion,  but 
also  modeling  it.  In  macroscopic  systems  con¬ 
taining  a  large  number  of  particles,  Fick’s  law  of 
diffusion  describes  the  time  evolution  of  particles 
moving  from  regions  of  high  to  low  concentration, 
resulting  in  a  smooth  (i.e.  fluctuation-free)  tran¬ 
sition  toward  equilibrium.  By  contrast,  polar¬ 
ization  in  few-spin  ensembles  can  instantaneously 
be  much  larger  than  the  Boltzmann  polarization, 
and  even  flow  “against”  the  polarization  gradi¬ 
ent  [12].  Our  measurements  are  corroborated 
with  a  model  which  numerically  simulates  this 
inherent  randomness  of  few-spin,  statistically- 
polarized  ensembles. 

Spin  centers  in  diamond  have  attracted  much 
attention  because  of  their  long  spin-lattice  relax¬ 
ation  time,  T\  [13].  In  order  to  study  flip-flop- 
dominated  spin  dynamics,  we  have  fabricated  a 
nano-scale  channel  of  implanted  PI  centers  in  di¬ 
amond  at  a  density  of  6  ppm.  The  flip-flop  fluctu¬ 
ations  are  thereby  constrained  to  this  ‘spin  wire.’ 
The  density  of  PI  centers  results  in  strong  dipo¬ 
lar  coupling  between  spins,  and  a  spin-spin  relax¬ 
ation  time,  X2  =  6.7  /is,  which  is  many  orders 
of  magnitude  shorter  than  Ti=  1.3  s.  Since  the 
flip-flop  time  Tff  is  equal  to  30  X2  [14],  the  spins 
in  the  wire  undergo  ~  104  flip-flops  for  each  T\ 
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relaxation  event. 

We  use  magnetic  resonance  force  microscopy 
(MRFM)  to  probe  the  spins  within  the  wire  (Fig. 
la).  The  force  detector  is  an  IBM-style  ultra-soft 
cantilever  [15]  with  a  SmCos  magnetic  particle 
glued  onto  the  tip.  To  resonantly  detect  spins,  we 
implement  the  iOSCAR  timing  protocol  [16].  A 
super-conducting  niobium  coil  is  used  to  drive  the 
spin  resonance  at  a  frequency  of  ^  =  2.18  GHz. 
This  frequency  defines  a  ‘resonant  slice’  where  the 
total  magnetic  field  (tip  field  plus  external  field) 
experienced  by  sample  spins  is  equal  to  ^  =  780 
G.  By  driving  the  cantilever  at  its  self-determined 
resonant  frequency  to  an  amplitude  of  150  nm, 
the  resonant  slice  sweeps  out  a  volume  of  the  same 
width,  hereafter  termed  the  detection  volume. 

We  scan  the  detection  volume  into  the  spin  wire 
and  measure  the  force  exerted  on  the  cantilever 
by  the  selected  spins.  From  this  measurement  we 
can  extract  two  quantities:  spin  signal  (Fig.  If, 
top),  and  correlation  time  (Fig.  If,  bottom).  The 
spin  signal  is  obtained  from  the  variance  in  the 
time  record  of  the  force  signal  and  is  directly  re¬ 
lated  to  the  number  of  measured  spins  [17].  The 
force  correlation  time,  rm,  describes  the  charac¬ 
teristic  time  for  the  net  moment  of  the  detected 
spins  to  decorrelate.  Once  the  detection  volume 
enters  the  channel,  the  spin  signal  grows  with  the 
number  of  detected  spins,  eventually  reaching  a 
plateau  when  completely  within  the  wire.  The 
correlation  time  shows  more  complex  behavior  as 
a  function  of  walk-in  distance.  These  results  can 
be  understood  within  the  framework  of  flip-flop- 
mediated  spin  diffusion:  ensemble  correlation  is 
lost  as  spin  diffuses  into  or  out  of  the  detection 
volume.  Just  inside  a  displacement  of  0  nm  in 
Fig.  lb,  the  force  signal  fluctuates  as  spins  in¬ 
side  the  volume  easily  and  rapidly  flip-flop  with 
outside  neighbors,  resulting  in  a  relatively  short 
correlation  time.  With  increasing  overlap  of  the 
detection  volume  and  the  spin  wire,  spin  must 
diffuse  further  to  exit  the  detection  volume  and 
change  the  overall  state  of  the  ensemble,  thus  in¬ 
creasing  rm.  The  correlation  time  peaks  when 
approximately  66  %  of  the  detection  volume  has 
entered  the  channel.  This  is  attributed  to  the 
force  sensitivity  profile  of  the  detection  volume: 
the  middle  of  the  volume  is  most  sensitive  to  spins 


(indicated  by  red  in  Fig  1),  while  the  edges  are 
least  sensitive  (blue;  see  Sec.  V).  As  the  detec¬ 
tion  volume  approaches  100  %  overlap  (at  150 
nm),  spins  are  able  to  diffuse  out  either  side  of 
the  most  sensitive  (red)  region,  reducing  the  cor¬ 
relation  time  by  roughly  a  factor  of  2.  Scan¬ 
ning  deeper  into  the  channel  results  in  no  fur¬ 
ther  change,  since  the  measurement  geometry  be¬ 
comes  translationally  invariant. 

A  global  measurement  of  the  spin  polarization, 
encompassing  the  entire  spin  wire,  would  not  un¬ 
cover  this  spin  transport  process.  The  correlation 
time  rm  plotted  in  Fig.  If  does  not  indicate  the 
time  scale  over  which  global  sample  magnetiza¬ 
tion  decays  away,  namely  T% |  rather,  the  correla¬ 
tion  time  is  dominated  by  spin  transport  out  of 
our  detection  volume.  As  such,  our  measurement 
is  similar  to  the  determination  of  spin-lifetime 
by  electrically-detected  three-terminal  Hanle  or 
optically-detected  “local  Hanle”  [18],  both  of 
which  can  be  obscured  by  spins  diffusing  away 
from  the  detection  volume  (electrical  contact  or 
optical  probe  spot).  However,  in  light  of  the  long 
T\  in  our  system,  the  flip-flop  mechanism  pro¬ 
vides  an  efficient  means  of  pure  spin  transport 
through  the  spin  wire.  In  the  limit  of  measuring 
a  single  spin,  Ti-type  relaxation  would  be  indis¬ 
tinguishable  from  flip-flop  induced  polarization 
changes.  We  can  distinguish  these  effects  because 
we  can  systematically  vary  the  size  of  the  over¬ 
lap  of  the  detection  volume  with  the  spin  wire,  as 
well  as  the  number  of  detected  spins. 

Since  the  polarization  gradient  vanishes  at 
thermal  equilibrium,  application  of  a  conven¬ 
tional  model  of  diffusion  driven  by  the  polariza¬ 
tion  gradient  would  lead  one  to  conclude  that 
rm  is  independent  of  detection  volume  position 
within  the  wire.  Such  a  model  neglects  the  spin 
fluctuations  which  dominate  for  small  ensembles. 
Therefore,  in  order  to  fit  the  measured  data  in 
Fig.  If,  we  use  a  Monte  Carlo  simulation  which 
simulates  the  flip-flops  between  individual  spins 
as  a  Markov  process.  From  this  fit  we  can  ex¬ 
tract  the  flip-flop  time  Tjj,  which  we  find  to  be 
0.21  ms.  Using  the  relationship  T2  =  =  6.7 

fi s,  we  find  good  agreement  with  previous  mea¬ 
surements  of  T2  (10.4  ±  7.1  /jls)  [19].  The  mea¬ 
sured  rm  can  be  viewed  as  the  time  needed  to 
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FIG.  1:  MRFM  spin  wire  setup,  scan  zoom-ins,  and  spin  wire  measurements,  a,  MRFM  setup 
for  measuring  the  200  nm  x  250  nm  x  4  /im  spin  wire  of  PI  centers  in  diamond.  The  coil  resonates  spins 
along  a  contour  of  the  magnetic  particle’s  field  called  the  ‘resonant  slice,’  and  as  the  cantilever  oscillates,  this 
slice  sweeps  out  a  spin  detection  volume.  Spins  toward  the  center  of  the  detection  volume  contribute  to  the 
spin  signal  more  strongly  (red)  than  spins  near  the  edge  (blue),  b-e,  Zoom-ins  showing  the  detection  volume 
scanned  into  the  stripe,  and  spins  diffusing  out.  f,  The  spin  signal  as  a  function  of  walk-in  distance  increases 
until  the  volume  is  fully  inside  the  spin  wire,  g,  The  correlation  time  of  the  measured  spins.  At  first  spins 
are  able  to  quickly  diffuse  out  of  the  volume,  but  as  the  measurement  volume  is  walked  in  further,  spins  must 
diffuse  a  longer  distance  to  exit  the  volume,  causing  the  correlation  time  to  increase.  Once  the  measurement 
volume  is  fully  inside  the  stripe,  spins  can  diffuse  out  both  sides,  decreasing  the  correlation  time  to  roughly 
half  the  peak  value. 
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FIG.  2:  Hyperfine  spectrum  measurement  for  nanowire,  a,  Energy  level  diagrams  for  three  spin 
detection  volumes  resulting  from  the  hyperfine  splitting.  The  three  nuclear  spin  populations  (raj  =  -1,  0, 
+1)  are  simultaneously  on  resonance,  but  due  to  the  magnetic  field  gradient,  this  occurs  at  three  distinct 
locations  in  the  sample,  b,  The  three  volumes  are  slowly  scanned  into  the  spin  wire,  one  at  a  time,  c,  The  spin 
signal  increases  in  a  step-like  fashion  as  each  volume  enters  the  wire  with  the  solid  line  showing  the  weighted 
smoothing  spile  used  for  signal  deconvolution  used  to  obtain  the  force  spectrum  shown  in  d.  d,  A  nanovolume 
EPR  spectrum  can  be  obtained  from  deconvolving  the  spin  signal  with  the  three  hyperfine  peaks  of  average 
FWHM  of  7.6  Gauss  resolved. 


transport  spin  out  of  the  detection  volume  via 
random- walk  diffusion  processes,  with  each  step 
taking  an  average  time  of  Tff.  From  the  mea¬ 
sured  flip-flop  rate,  we  can  extract  a  diffusion 
constant  of  D  =  ^  4.5  x  10-9  where 

Tff  s 

a  =  9.82  nm  is  the  average  nearest-neighbor  sep¬ 
aration  for  spins  in  the  6  ppm  wire.  The  diffusion 
length  for  a  pure  spin  current  can  then  be  deter¬ 
mined  by  L  =  y/DTi,  for  which  we  find  721  nm. 
This  is  competitive  with  metallic  spin  transport 
devices  at  room  temperature  [20]. 

Suppression  of  diffusion  can  occur  when,  due 
to  a  magnetic  field  gradient,  neighboring  spin 


sites  experience  different  Zeeman  splitting  which 
would  cause  flip-flops  to  violate  energy  conser¬ 
vation  [21,  22].  This  mismatch  in  Zeeman  split¬ 
ting  can  be  compensated  if  there  is  substantial 
variation  of  the  spin-split  levels  throughout  the 
sample.  This  variation  is  characterized  by  the 
ESR  linewidth.  We  measure  an  inhomogeneously 
broadened  linewidth  of  7.6  Gauss  as  discussed 
below  and  in  Sec.  IV.  This  linewidth  is  signif¬ 
icantly  broader  than  those  which  have  been  ob¬ 
served  in  other  MRFM  spectroscopy  experiments, 
and  therefore  we  do  not  observe  suppression  of 
diffusion  in  our  experiment. 
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Next,  we  repeat  the  experiment  with  a  vari¬ 
ation  of  the  iOSCAR  protocol  in  order  to  si¬ 
multaneously  improve  spatial  resolution  and  pre¬ 
serve  maximum  force  sensitivity.  A  conventional 
iOSCAR  measurement  utilizes  the  entire  detec¬ 
tion  volume  swept  out  by  the  cantilever  oscilla¬ 
tion  (red  and  blue  regions,  Fig.  lb-e),  as  this 
couples  to  the  largest  number  of  spins  and  cre¬ 
ates  the  largest  force  signal.  By  truncating  the 
detection  volume  to  the  most-sensitive  (red)  por¬ 
tion,  we  couple  only  to  spins  in  this  smaller  region 
while  maintaining  high  force  sensitivity.  In  con¬ 
trast,  reducing  the  cantilever  oscillation  ampli¬ 
tude,  while  similarly  improving  resolution,  sim¬ 
ply  compresses  the  sensitivity  profile  by  shrinking 
the  red  region  and  retaining  the  reduced  sensitiv¬ 
ity  near  the  edges.  Utilizing  this  technique  (see 
Sec.  V),  which  we  call  partially-interrupted  OS¬ 
CAR  (piOSCAR),  we  again  scan  the  detection 
volume  into  the  spin  wire.  The  improved  reso¬ 
lution  enables  us  to  resolve  a  staircase  structure 
with  three  steps  (Fig.  2c),  corresponding  to  hy- 
perfine  splitting  as  discussed  below. 

The  PI  center  is  a  substitutional  nitrogen  de¬ 
fect  with  an  unpaired  spin-^  electron  and  spin- 
1  nucleus.  Due  to  hyperfine  interaction,  each 
Zeeman-split  electron  energy  level  splits  into 
three  states  (Fig.  2a)  that  correspond  to  the 
three  states  for  the  nuclear  spin,  mj  =  —1,  0,  1. 
The  hyperfine  splitting  for  our  diamond  crystal 
orientation  with  respect  to  the  external  field  is 
33  G  [23].  In  the  approximately  dipolar  field  of 
our  probe  magnet,  this  splitting  manifests  itself 
as  three  non-overlapping  spin  detection  volumes, 
with  each  volume  defined  by  the  contour  of  mag¬ 
netic  field  which  satisfies  the  resonance  condition 
for  one  hyperfine  transition  (Fig.  2b).  We  ob¬ 
serve  a  step-like  increase  in  the  spin  signal  as 
each  volume  enters  the  wire  and  the  cantilever 
becomes  coupled  to  an  additional  hyperfine  tran¬ 
sition.  The  broadening  of  the  steps’  rising  edges 
is  caused  by  the  width  of  each  volume  (b  =  15 
nm).  The  spacing  between  step  edges  (s  =  25 
nm),  along  with  the  hyperfine  splitting  constant 
(33  G)  provides  a  means  of  measuring  the  MRFM 
probe  field  gradient,  G  =  1.3  G/nm.  The  stair¬ 
case  structure  is  a  convolution  of  the  implanted 
channel  profile  (essentially  a  step  function),  our 


force  sensitivity  profile,  and  the  intrinsic  ESR 
spectrum  of  spins  in  the  sample.  We  can  there¬ 
fore  obtain  the  ESR  spectrum  through  the  ap¬ 
propriate  deconvolution  steps  (see  Sec.  IV).  This 
spectrum  is  shown  in  Fig.  2d. 

Our  measurements  demonstrate  that  both  spin 
transport  and  spectroscopic  information  can  be 
obtained  by  observing  spin  fluctuations  about  the 
equilibrium  in  a  nano-scale  diamond  spin  wire. 
We  are  able  to  distinguish  spin  transport  from 
Ti-type  spin  relaxation.  We  find  a  flip-flop  time 
Tff  of  0.21  ms  and  determine  a  spin  diffusion 
constant  of  D  ^  4.5  x  10-9  for  6  ppm  PI 
centers  in  diamond.  This  regime  —  small  spin 
ensembles  which  experience  substantial  spin-spin 
coupling  —  could  prove  useful  for  the  develop¬ 
ment  of  spintronic  technologies,  especially  as 
technology  trends  towards  ever-smaller  devices. 
Additionally,  the  capability  to  spatially  resolve 
nuclear  spin  states  might  allow  for  the  selective 
study  of  nuclear  spin  dynamics  via  electron 
spin  detection  in  a  manner  reminiscent  of  bulk 
electron-nuclear  double  resonance  (END OR), 
but  at  the  nano-scale.  In  all,  our  measurements 
elucidate  the  mechanisms  of  spin  transport  in 
small  ensembles,  and  are  a  valuable  tool  for  the 
study  of  nano-scale  spin  phenomena  and  devices. 

Methods  The  ultra-soft  silicon  cantilever  has 
a  spring  constant  of  100  fi N/m,  resonance  fre¬ 
quency  of  about  3  kHz,  and  approximate  dimen¬ 
sions  of  90  microns  in  length,  1  micron  width, 
and  100  nm  thickness.  Displacement  of  the  can¬ 
tilever  is  measured  through  a  1550  nm  laser  inter¬ 
ferometer.  The  MRFM  measurements  were  taken 
at  temperatures  of  4.2  K.  The  measured  moment 
of  the  SmCos  magnetic  particle  was  3.6  x  10-12 
J/T.  The  niobium  resonator  coil  had  about  2.5 
turns  and  a  300  micron  diameter,  and  the  reso¬ 
nance  bowl  had  a  diameter  of  about  3  microns. 
The  diamond  sample  studied  in  this  experiment 
had  a  background  nitrogen  concentration  of  0.3 
ppm  (5.27  x  1016  cm-3).  To  create  the  high  spin- 
density  wire,  the  sample  was  first  prepared  using 
electron  beam  patterning  and  then  exposed  to  ni¬ 
trogen  ion  implantation  at  a  variety  of  energies  to 
create  a  uniform  spin  density.  Finally  the  sample 
was  annealed  to  yield  an  approximately  uniform, 
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channel  (6  ppm  or  1.06  x  1018  cm-3)  with  width, 
depth,  and  length  of  200  nm,  250  nm,  and  4  mi¬ 
crons,  respectively.  The  sample  was  purchased 
commercially  from  Element  Six,  grown  via  chem¬ 
ical  vapor  deposition,  and  the  nitrogen  ion  im¬ 
plantation  was  performed  by  Leonard  Kroko  Inc. 


III.  MICROSCOPIC  MONTE  CARLO 

SIMULATION  OF  SPIN  DIFFUSION  IN 
QUASI  ID  SYSTEMS 

To  fit  our  data,  we  developed  a  Monte  Carlo 
simulation  of  spins  flip-flopping  according  to  the 
random-walk  dynamics  as  originally  introduced 
by  Bloembergen,  and  discussed  by  Abragam  [14]. 
In  the  experiment,  we  are  only  sensitive  to 
changes  in  magnetization  along  the  x-direction, 
defining  x  to  be  along  the  wire,  which  led  us  to 
implement  a  one-dimensional  simulation  for  com¬ 
putational  efficiency.  In  our  simulation  we  break 
up  the  spin  wire  into  an  array  of  spins,  then  let 
the  spins  randomly  walk.  We  make  the  approx¬ 
imation  that  all  spins  can  flip  with  neighbors 
that  are  located  at  an  arbitrary  0  and  0,  refer¬ 
ring  to  spherical  coordinates,  and  take  the  pro¬ 
jection  along  the  x-axis  for  each  flip-flop  event; 
thus  if  an  up  spin  flip-flops  to  a  nearest  neighbor 
(6  ppm  o  9.82  nm  nearest  neighbor  separation) 
at  an  angle  of  6  =  45°,  0  =  0,  then  we  are  only 
interested  in  the  x-component  of  r sin# cos 0  = 
9.82  sin  (45°)  cos  (0)  =  9.82/a/2  nm,  which  we  use 
to  calculate  the  change  in  the  spatial  distribu¬ 
tion  of  spin  (magnetization).  This  approxima¬ 
tion  is  valid  in  the  limit  that  the  discrete  0  and 
0  imposed  by  the  crystal  lattice  become  negli¬ 
gible;  since  the  lattice  constant  for  diamond  is 
3.57  A,  adjacent  lattice  sites  differ  by  about  2  de¬ 
grees  (from  a  distance  away  of  9.82  nm),  which 
we  consider  a  discretization  that  is  small  enough 
to  ignore.  A  further  consideration  with  this  ap¬ 
proximation  is  that  spins  on  the  edges  of  the  wire 
have  less  neighbors,  and  thus  cannot  flip-flop  at 
any  arbitrary  0  or  0.  As  a  consequence  of  having 
less  neighbors,  these  edge  spins  flip-flop  less  fre¬ 
quently,  and  since  there  are  significantly  less  of 
them  than  interior  spins,  we  also  consider  their 
effects  small  enough  to  ignore. 


The  random  walk  dynamics  in  our  model  re¬ 
quire  a  step  size,  and  a  rate  at  which  to  take  the 
steps.  The  step  size  for  our  system  is  the  average 
nearest-neighbor  separation,  which  has  large  vari¬ 
ation  given  the  fact  that  the  nitrogen  atoms  end 
up  on  random  lattice  sites  during  implantation. 
Despite  this  variation  in  nearest-neighbor  sepa¬ 
ration,  which  causes  more  closely-packed  spins  to 
flip-flop  more  rapidly,  and  more  isolated  spins  to 
flip-flop  less  rapidly,  the  average  flip-flop  rate  for 
spins  is  given  by  W  =  =  3.2  kHz.  This 

gives  the  average  flip-flop  time  Tff  =  for  the 
time-random  flip-flop  transitions  modeled  in  the 
simulation.  This  treatment  of  the  flip-flop  time 
avoids  the  computationally  prohibitive  alterna¬ 
tive  of  calculating  the  flip-flop  rate  between  all 
pairs  of  spins. 

Conceptually,  we  represent  a  flip-flop  event  by 
having  a  single  spin  hop,  instead  of  having  two 
neighboring  spins  undergo  simultaneous  transi¬ 
tion  flips,  which  is  a  way  of  increasing  compu¬ 
tational  efficiency  by  using  half  the  number  of 
spins  to  accurately  represent  the  magnetization 
record  of  all  the  spins.  This  can  be  thought  of  as 
shifting  the  value  of  magnetization  for  all  spins 
up  by  a  Bohr  magneton  —  setting  down  spins  to 
having  zero  moment,  and  up  spins  to  2  /ib  —  be¬ 
cause  when  an  actual  flip-flop  occurs,  there  is  the 
transport  of  2  fiB  along  the  direction  of  the  vec¬ 
tor  connecting  the  two  spins.  Whereas  our  mea¬ 
surement  fluctuates  about  zero  net  moment  with 
roughly  y  spins  up  and  y  spins  down  (where 
N  is  the  number  of  spins  that  we  measure),  our 
simulation  fluctuates  about  an  average  moment 
of  y  x  (2 /ib)  =  At//#,  which  is  an  offset  that  we 
subtract  off  to  obtain  the  magnetization  record. 

To  simulate  the  dynamics,  we  generate  a 
Markov  matrix  that  contains  the  probabilities 
that  a  given  spin  transitions  from  a  particular 
spatial  bin  (represented  by  a  column  of  the  ma¬ 
trix)  to  a  particular  spatial  bin  (rows  of  the  ma¬ 
trix)  in  a  time  interval  At.  For  example,  a  spin 
in  the  middle  of  the  wire  can  flip-flop  a  distance 
of  9.82  nm  in  any  direction,  at  an  average  rate 
W,  corresponding  to  a  total  probability  of  hop¬ 
ping  in  a  time  At  of  P  =  W At.  To  eliminate 
the  chance  of  an  individual  spin  undergoing  two 
or  more  flip-flops  in  the  time  interval,  the  con- 
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FIG.  3:  Schematic  illustration  of  the  simulation.  As 
depicted,  the  detection  volume  has  walked  in  three 
steps,  and  measures  the  first  three  bins  of  the  spin 
array.  The  spins  randomly  walk  throughout  the  array 
resulting  in  time  dependent  total  force.  Autocorrela¬ 
tion  of  the  force  gives  the  spin  correlation  time  for 
the  given  location  of  the  detection  volume. 


dition  P  «  1  must  be  satisfied.  We  take  only 
the  x-component  of  the  hop,  which  must  be  dis¬ 
cretized  into  bins:  given  that  the  spins  can  hop 
in  any  direction  with  equal  probability,  the  prob¬ 
ability  distribution  for  the  x-component  is  a  flat 
distribution  out  to  the  nearest-neighbor  separa¬ 
tion.  Thus,  for  constructing  the  elements  (prob¬ 
abilities)  of  the  Markov  matrix  we  have  for  each 
element  one  of  three  cases:  1.)  the  value  of  the 
element  is  0  because  it  belongs  to  a  spatial  bin 
that  is  too  far  away  (>  9.82  nm)  for  the  spin  to 
hop,  2.)  the  element  is  equal  to  p  =  where 
Nib  =  (  2x9.88 nm — \  -g  num|3er  of  local 

bins  that  are  ’within  reach’,  3.)  the  element  is 
on  the  diagonal  and  is  thus  1  minus  the  sum  of 
the  other  values  in  the  column.  The  values  of  the 
diagonals  are  chosen  such  that  the  columns  are 
properly  normalized,  where  the  total  probability 
of  the  spin  to  go  somewhere  (including  staying  in 
its  current  bin)  is  unity. 

The  magnetization  array  divides  the  wire  into 
a  number  of  bins  of  finite  size  Ax,  such  that  there 
are  only  a  few  net  spins  in  each  bin.  As  an  ap¬ 
proximation,  we  take  the  center  of  the  bin  to  be 
the  location  of  all  the  spins  in  the  region  that 
the  bin  represents,  which  becomes  more  accurate 


as  the  bin  size  gets  smaller.  We  let  the  spins 
in  our  array  propagate  in  time  according  to  the 
described  dynamics,  generating  a  new  magneti¬ 
zation  record  for  each  small  time  step  At.  To 
simulate  our  experimental  measurement  of  walk¬ 
ing  into  the  wire,  we  multiply  the  magnetization 
array  by  our  spin  sensitivity  profile  (as  described 

in  Sec.  V,  with  units  of  m^ment ) :  in  the  first 
step  we  multiply  the  first  bin  in  our  magneti¬ 
zation  array  by  the  leading  edge  of  the  sensi¬ 
tivity  profile,  followed  by  multiplying  the  first 
two  bins  by  the  corresponding  sensitivity,  and 
so  on  and  so  forth,  resulting  in  a  force  record 

( moment  X  momen^  =  f°rce)  for  each  step  into 
the  wire.  To  extract  the  correlation  time,  we  take 
the  autocorrelation  of  the  force  record  and  fit  a 
decaying  exponential;  the  exponential’s  time  con¬ 
stant  is  the  correlation  time.  Fitting  a  decaying 
exponential  to  the  autocorrelation  assumes  that 
the  flip-flops  can  be  modeled  as  a  Poisson  random 
process  where  the  probability  P  for  a  flip-flop  to 
occur  is  the  same  for  every  time  step  At.  Typi¬ 
cal  parameters  for  the  simulation  used  to  fit  the 
data  are:  Number  of  At  time  steps  r\j  107,  flip-flop 
probability  P  ~  0.02,  and  Ax  ~  3  nm. 


IV.  SPECTROSCOPIC 
MEASUREMENTS  ON  NANOSCALE 
VOLUME 

The  data  presented  in  Fig.  2c  shows  the  ES- 
RFM  force  recorded  as  the  cantilever  with  a 
probe  magnet  is  brought  closer  to  the  end  of 
the  ion  implanted  spin  nanowire.  The  distinctive 
staircase-shaped  pattern  of  increasing  force  sig¬ 
nal,  observed  in  the  data,  is  consistent  with  the 
three  individual  sensitive  volumes,  defined  by  the 
splitting  of  the  ESR  line  by  the  hyperfine  interac¬ 
tion  given  by  Eq.  1,  entering,  consecutively,  the 
end  of  the  nanowire.  The  obtained  data,  com¬ 
bined  with  our  knowledge  of  the  spatial  distribu¬ 
tion  of  the  spins  in  the  nanowire  and  the  spatial 
width  of  the  rf  interrupt  (see  Sec.  V)  allows  us  to 
extract  spectroscopic  information  about  the  spins 
in  the  nanowire. 

The  Hamiltonian  for  an  electron  on  a  PI  cen¬ 
ter  has  four  terms:  electron  Zeeman,  hyperfine 
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a) 


A fhf 
G 


LI  I 


—  Force  spectrum 


b) 


Point  Spread  Function 


Awint  x  G  ,  Awint  +  -g 

Interrupt  Point  Spread  Function 


Spin  Distribution  Observed  Force 


FIG.  4:  Schematic  representation  of  a  two-step  for¬ 
mation  of  the  force  signal  shown  in  Fig.  2c.  a)  A  hy- 
perfine  split  spectrum  of  an  electron  on  a  PI  center 
results  in  spatial  force  spectrum  consisting  of  three 
sensitive  slices  of  spatial  width  of  A  A  /G,  where  A  A  is 
the  ESR  linewidth  and  G  is  the  magnetic  field  gradi¬ 
ent,  which  are  spatially  separated  by  A Hhf  /G,  where 
A Hhf  is  the  hyperfine  splitting  of  the  spectrum.  The 
force  spectrum  convolved  with  the  shape  of  the  in¬ 
terrupt,  incorporating  cantilever  oscillation  and  mi¬ 
crowave  manipulation,  as  described  in  Sec.  V,  forms 
the  Point  Spread  Function  (PSF)  of  the  ESRMRF.  b) 
The  convolution  of  the  PSF  with  the  step  function¬ 
like  spatial  distribution  of  the  spins  in  the  nanowire 
results  in  the  observed,  staircase  shaped  force  pattern 
shown  in  Fig.  2c. 


interaction,  quadrupole,  and  nuclear  Zeeman, 

T-L  =  gePBS  ■  H  +  a£X  +  PIl  ~  9rUhJ-  H 

Electron  Hyperfine  Quadrupole  Nuclear 
Zeeman  Zeeman 

(i) 

where  ge  and  gn  are  the  electron  and  nitrogen 
nucleusg-factors,  /ab  and  gn  are  the  Bohr  mag¬ 
neton  and  nuclear  magneton  respectively,  S  and 
/  are  the  spin  of  the  electron  and  nucleus,  and 
H  is  an  externally  applied  magnetic  field.  The 
hyperfine  term  of  the  Hamiltonian  is  responsible 
for  the  splitting  of  the  ESR  spectrum  into  three 
lines  separated  by  the  hyperfine  splitting  A  Hhf 
=  33  Gauss. 

The  general  principle  of  the  formation  of  the 
observed  force  pattern  is  shown  schematically  in 
Fig. 4.  The  force  pattern  formation  occurs  in  two 
steps.  A  hyperfine  split  spectrum  of  an  electron 
on  a  PI  center  results  in  spatial  force  spectrum 


FIG.  5:  a)  ESRFM  point  spread  function  (PSF)  re¬ 
covered  via  deconvolution  of  the  force  signal  shown  in 
Fig.  2c  with  the  step  function-like  spatial  distribution 
of  the  spins  in  the  nanowire,  b)  Spatial  representa¬ 
tion  of  the  force  spectrum  of  electron  spins  on  PI 
center  in  diamond  recovered  via  deconvolution  of  the 
PSF  with  the  shape  of  a  15  nm  interrupt. 


(see  Fig. 4a)  consisting  of  three  sensitive  slices  of 
spatial  width  of  AA/G,  where  A  A  is  the  ESR 
linewidth  and  G  is  the  magnetic  field  gradient, 
which  are  spatially  separated  by  AH^f/G,  where 
A  Hhf  is  the  hyperfine  splitting  of  the  spectrum. 
The  force  spectrum  convolved  with  the  shape 
of  the  interrupt,  incorporating  cantilever  oscilla¬ 
tion  and  microwave  manipulation,  as  described  in 
Sec.  V,  forms  the  Point  Spread  Function  (PSF) 
of  the  ESRMRF.  Subsequently,  the  convolution 
of  the  PSF  with  the  step  function-like  spatial  dis¬ 
tribution  of  the  spins  in  the  nanowire  results,  as 
illustrated  in  Fig. 4b,  in  the  observed,  staircase 
shaped  force  pattern  shown  in  Fig.  2c. 

Based  on  the  described  model,  we  have  ex¬ 
tracted  the  shape  of  the  electron  spin  spectrum  in 
the  implanetd  nanowire  using  the  known  shape  of 
the  interrupt  and  the  good  knowledge  about  the 
spatial  distribution  of  spins  in  the  nanowire  with 
the  result  shown  in  Fig.  5b.  This  was  achieved 
via  two-step  deconvolution  process.  In  the  first 
step,  the  experimental  data  approximated  by  a 
weighted  smoothing  spine  fit,  shown  in  Fig.  2c,  is 
deconvolved  with  a  step  step  function,  approxi¬ 
mating  the  spatial  spin  distribution,  resulting  in 
the  PSF  shown  in  Fig.  5a.  The  deconvolution  was 
done  numerically,  using  a  modified  Van-Cittert  it¬ 
erative  procedure  of  deconvolution.  During  the 
second  step,  the  obtained  PSF  is  deconvolved 
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with  the  interrupt  shape  approximated  by  a  15 
nm  wide  square  pulse  resulting  in  the  force  spec¬ 
trum  shown  in  Fig.  5b.  This  deconvolution  was 
done  using  iterative  Bayesian-based  method  [24]. 
The  resulting  spectral  linewidth  of  7.6  Gauss,  is 
attributed  to  the  inhomogenious  line  broadening 
and  an  therefore  contributes  to  the  explanation 
of  the  absence  of  suppression  of  spin  diffusion  in 
our  experiment  even  in  the  presence  of  the  strong 
magnetic  field  gradient. 


iOSCAR  sensitivity  profile  we  can  then  show  how 
implementing  partial  interrupts  in  pi-OSCAR  re¬ 
sults  in  greater  spin  sensitivity  per  unit  volume 
for  small  detection  volumes. 

We  treat  the  cantilever  in  our  experiment  as  a 
simple  damped,  driven  harmonic  oscillator  with 
spring  constant  k ,  viscous  damping  coefficient  c, 
effective  mass  m,  and  driving  force  F(t).  From 
Newton’s  second  law  we  have 


V.  PARTIALLY  INTERRUPTED  IOSCAR 
PROTOCOL 


mx  =  F(t)  —  kx  —  cx  (2) 

which  can  be  rewritten  in  a  more  conventional 
form: 


We  implemented  a  slightly  modified  version  of 
the  iOSCAR  protocol  which  improves  spatial  res¬ 
olution,  which  we  call  partially-interrupted  oscil¬ 
lating  cantilever-driven  adiabatic  reversals,  or  pi- 
OSCAR.  The  iOSCAR  protocol  involves  periodi¬ 
cally  interrupting  the  microwave  power  for  a  half 
cantilever  cycle,  which  causes  the  orientation  of 
spins  in  the  measurement  volume  to  change  by 
180  degrees  with  respect  to  the  cantilever  mo¬ 
tion.  This  in  turn  causes  the  sign  of  the  force  on 
the  cantilever  to  be  periodically  reversed,  which 
can  then  be  demodulated  via  a  lock-in  amplifier 
to  obtain  the  spin  signal.  The  use  of  a  lock-in 
allows  for  increased  spin  sensitivity.  The  modifi¬ 
cation  we  implement  in  pi-OSCAR  is  to  interrupt 
the  microwave  power  for  less  than  a  half  can¬ 
tilever  cycle,  about  the  center  of  the  cantilever’s 
motion.  Since  only  the  spins  that  are  periodi¬ 
cally  interrupted  contribute  to  the  spin  signal  due 
to  the  lock-in  detection,  by  decreasing  the  width 
of  the  interrupt  it  correspondingly  decreases  the 
measurement  volume  width.  Using  this  modifi¬ 
cation,  we  get  larger  signals  from  the  same  mea¬ 
surement  volume  compared  to  iOSCAR,  enabling 
us  to  probe  down  to  smaller  measurement  vol¬ 
umes  before  we  become  limited  by  noise.  The 
larger  signal  is  due  to  the  fact  that  the  center  of 
the  spin  detection  volume  has  a  larger  force  sen¬ 
sitivity  per  spin  than  the  edges  of  the  volume. 
The  following  derivation  solves  for  the  force  sen¬ 
sitivity  of  the  detection  volume  as  a  function  of 
lateral  position  within  the  volume  for  the  stan¬ 
dard  iOSCAR  protocol  measurement.  From  this 


x  +  2Toj0x  +  ojqX 


m 

m 


(3) 


where  the  damping  ratio  is  T  = 
undamped  angular  frequency  is 


c 

2  \Jmk 


and  the 


w0  =  2tt/  =  J—  (4) 

V  m 

Driving  the  cantilever  results  in  a  steady  state 
oscillation  such  that  the  energy  dissipated  each 
cycle  is  supplied  by  F(t)  (in  accordance  with  its 
quality  factor).  Thus  the  position  of  the  can¬ 
tilever’s  magnetic  tip  is  given  by 


xc(t)  =  xpk  sin (uct)  (5) 

where  ujc  =  27r/c  is  the  steady  state  cantilever  fre¬ 
quency  and  xpk  is  the  cantilever  peak  oscillation 
amplitude.  Now,  as  the  cantilever  oscillates,  the 
resonance  bowl  sweeps  past  spins  within  the  spin 
detection  volume  and  causes  them  to  adiabati- 
cally  invert.  This  results  in  the  spin’s  magnetic 
moment  rh  to  vary  as  a  pulse  wave  in  time: 


m(t)  =  fJ,Bfpulse(t)  (6) 

here,  hb  is  a  Bohr  magneton  and  fpuise(t)  is  a 
pulse  wave  modultion  signal  with  a  pulse  period 
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FIG.  6:  Pulse  function  of  two  spins’  moments  at  dif¬ 
ferent  locations  in  the  detection  volume. 


(width)  Tpuise  and  unit  amplitude: 


only  the  Fourier  series  component  of  Fspin(t)  at 
the  cantilever  frequency  ujc.  This  gives: 


—XpkUJ^  sm(u:ct)  +  2YxpkuJouoc  cos  (c uct)  + 

k  .  .  F(t)  GnB ,  ,  x  (9) 

— xpksm (uct)  = - 1 - bi  sm(uct) 

m  mm 

where  b\  is  the  dimensionless  amplitude  of  the 
Fourier  component  of  interest  (defined  below). 
Note  that  the  3rd  term  on  the  1.  h.  s.  can  be 
grouped  with  the  2nd  term  on  the  r.  h.  s.: 


/puZse(^) 


1 


for  T' 


Tpuls 


<  t  <  T'  + 


Tpuls 


—  1  otherwise 


(7) 

where  T'c  =  Tc(\  +  m)  and  m  =  0, 1,2,3...  and 
Tc  =  -j-  is  the  cantilever  period.  For  a  spin  in 
the  center  of  the  detection  volume  this  results  in 
a  50%  duty  cycle  pulse  wave  (i.e.  a  square  wave 
with  Tpuise  =  F-)  modulation  of  the  spin’s  mo¬ 
ment  fh  with  a  period  equal  to  that  of  the  can¬ 
tilever,  and  for  spins  not  in  the  center,  the  pulse  is 
greater  or  less  than  a  50%  duty  cycle  (see  Fig.  6). 
Because  it  is  the  cantilever’s  magnetic  tip  that 
cause  the  spins  to  invert,  the  pulse  wave  modu¬ 
lation  is  necessarily  phase-locked  and  frequency- 
locked  to  the  cantilever  motion.  This  modulation 
of  m  provides  a  small  additional  driving  force 
on  the  cantilever  Fspin(t )  =  | ^ (rh(t)  •  B)\  = 
G^b  f pulse  (t)  ,  where  B  is  the  magnetic  field  from 
the  cantilever  tip  and  G  is  appropriate  magnetic 
field  gradient,  which  we  take  to  be  uniform  within 
the  shallow  spin  wire.  Adding  this  additional 
force  to  eq.  9  and  plugging  in  the  steady  state 
solution  from  eq.  5  yields: 


—Xpkool  sin (cjct)  +  2Yxpk(jJoUJc  cos(u;ct)  + 


ujQXpk  sin (fjjct)  = 


m 


m 


Fspinjt )  (8) 

m 


but  since  Fspin(t)  is  a  pulse  wave  modulation  that 
is  frequency  and  phase  locked  to  the  oscillation 
of  the  cantilever,  and  because  the  cantilever  itself 
shows  sensitive  response  only  near  its  resonant 
frequency  (for  high  quality  factor),  we  can  keep 


—xpkuil  sm(wct)  +  2Txpkuj0ujc  cos(cjct)  + 

Gusbi  xpk  .  F(t)  (10) 

[k - — —  sm  (<jJct)  =  - 

xpk  m  m 

This  demonstrates  that  the  modulation  of  the 
spin  effectively  results  in  a  small  shift  5k  in  the 
cantilever’s  spring  constant: 


6k  =  (ii) 

Tpk 

From  the  undamped  angular  frequency  eq.  4,  we 
find  for  small  changes  in  the  spring  constant  5k: 


and  upon  substitution  back  into  eq  4  this  gives  a 
well  known  result  which  shows  that  small  changes 
in  the  spring  constant  give  a  linear  response  to 
changes  in  the  oscillator’s  frequency: 


5f  _  15k 

J~  2¥ 


(13) 


Now  we  need  to  find  the  Fourier  component  b\ 
for  a  given  spin  (not  necessarily  in  the  center  of 
the  volume).  First,  we  define  the  Fourier  series 
of  a  given  function  f(t)  with  period  T: 
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fit)  =  a0  +  J2 


,  2iTnt  v  ,  2iTnt  N 

cos(-— )  +  bn  sm(— — — ) 


n=l 

and  Fourier  coefficients: 


(14) 


T/2 

«o  =  ^  f  f(t)dt  (15) 

-T/2 


T/2 

2  /*  .  ,2itnt.  _  .  . 

«n  =  y  J  /(*)cos(-^-)dt  (16) 

-T/2 

T/2 

7  2  f  n  /  \  .  /27t nts  7  .  _ 

(17) 

-T/2 

where  /(£)  in  this  case  is  the  pulse  modulation  of 
the  spin  fpui8e(t )•  By  plugging  eq.  7  into  eq.  17 
and  calculating  a  straightforward  integration  we 
can  solve  for  &i  as  a  function  Tpuise : 

^1  impulse) 

T/2 

2  f  /27rF  . 

=  T  J  7(*)sm(  — )eft  (18) 

-T/2 

4  .  /  TpUise  \ 

=  -sm(— — tt) 

7 r  1  c 


FIG.  7:  First  harmonic  Forier  series  amplitude  &i,  as 
a  fuction  of  moment  pulse  width. 


gives  tfnp  =  |^arcsin(^).  However,  the  func¬ 
tion  arcsin(x)  has  a  restricted  range  —  |  <  x  <  | 
as  a  consequence  of  satisfying  the  vertical  line 
test  for  a  well  defined  (one-to-one  mapped)  func¬ 
tion.  As  a  consequence  of  this  we  need  to  pick  up 
the  ’other’  solution  (sometimes  referred  to  as  the 
geometric  reflection),  and  the  spacing  between 
these  solutions  gives  the  pulse  width.  Since  the 
other  solution  is  symmetric  about  |  ,  where  | 
corresponds  to  one  fourth  of  the  cantilever  pe¬ 
riod,  we  can  write  the  pulse  width  as  Tpuise(xs )  = 

2(x  -  i#  arcsin&))  =  Tail  -  1  arcsin&))- 


Tpulsei^  s) 


Tc  -  (  Xs  w 

—  arcsm( - )) 

2  7T  %pk 

f  •  /  XS  \\ 

—  arcsm( - )) 

7T  %pk 


(19) 


Fig.  7  shows  6i  as  a  function  of  Tpuise. 

Having  found  bi(Tpuise ),  we  now  need  to  re¬ 
late  the  pulse  wave  of  a  given  spin  to  its  lat¬ 
eral  location  xs  within  the  detection  volume.  Let 
xs  =  0  be  the  center  of  the  detection  volume  (i. 
e.  the  location  of  the  resonant  slice  when  the 
cantilever  is  at  its  equilibrium  position).  From 
eq.  5  the  range  of  xs  is:  —  xpk  <  xs  <  xpk- 
We  note  that  the  time  tfup  at  which  the  spin 
will  flip  as  a  function  of  its  location  xs  can  be 
found  by  equating  xs  =  xpk  sm(27rf ctfup)  which 


Having  solved  for  Tpuise  as  a  function  of  spin  lo¬ 
cation  xs,  we  can  now  combine  eqs.  18  and  19  to 
give  b\  (xs): 

7  .  4  .  7T  .  .  Xs  , 

bi(xs)  =  —  sm( - arcsm( - ))  (20) 

7T  2  Xpk 

Since  the  b\  coefficient  of  a  detected  spin  is  pro¬ 
portional  to  the  frequency  shift  of  the  cantilever 
due  to  that  spin  (eqs.  11  and  13)  we  have: 
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FIG.  8:  Force  sensitivity  profile  of  spin  detection  vol¬ 
ume  as  a  function  of  lateral  location. 


FIG.  9:  Comparison  of  iOSCAR  and  piOSCAR  rf 
pulse  signals  and  the  resultying  force  sensitivity  pro¬ 
file. 


5f(xs)  =  sin -  arcsin(  — ))  (21) 

7T  %pk  2  %pk 

We  see  that  6i(0)  =  ^  for  a  spin  in  the  center  of 
the  detection  volume,  which  gives  the  frequency 
shift: 


2  G^Bfc 

7 r  kxpk 


(22) 


Although  we  use  frequency  detection  in  our  mea¬ 
surements  and  therefore  Sf  is  a  useful  quantity,  it 
is  often  more  useful  to  report  the  corresponding 
forces  that  are  detected  from  the  spins  acting  on 
the  cantilever.  Therefore,  we  can  convert  back  to 
force  via  Hooke’s  law  (F  =  5k  •  xpk)  where  xpk 
is  used  to  yield  the  peak  force  and  5k  is  given  by 
eq.  11  which  gives  the  force  sensitivity  per  spin 
of  the  detection  volume  as  a  function  of  lateral 
position. 


F  ( xs )  =  GfiB 


4  .  ( 7T 

-sin - arcsm 

|_tt  V2 


^  1  aN 


(23) 

Having  found  the  force  sensitivity  of  the 
iOSCAR  protocol,  we  can  consider  the  modifi¬ 
cation  of  this  in  pi-OSCAR.  The  partial  inter¬ 
ruptions  in  the  pi-OSCAR  protocol  result  in  a 
truncated  force  sensitivity,  since  only  spins  that 
are  interrupted  contribute  to  force  detection  (see 
Fig.  9).  This  provides  an  oportunity  for  in- 


piOSCAR 

iOSCAR 

Decreasing 

Detection 

Volume 

A 

i 

A 

FIG.  10:  Comparison  of  iOSCAR  and  piOSCAR  for 
decreased  detection  volume  width. 


creasing  the  force  sensitivity  of  small  detection 
volumes.  Typically,  in  the  iOSCAR  protocol, 
smaller  volumes  are  probed  by  decreasing  the 
cantilever  oscillation  amplitude,  and  thus  shrink¬ 
ing  the  entire  force  sensitivity  profile  (Fig.  10). 
With  pi-OSCAR,  the  canntilever  oscillation  am¬ 
plitude  is  kept  at  a  large  fixed  value,  and  the 
volume  is  decreased  by  applying  partial  interrup¬ 
tions  to  the  microwave  resonator,  and  thus  trun¬ 
cating  the  sensitivity  profile  (Fig.  10).  This  re¬ 
sults  in  relatively  large  force  sensitivity  for  small 
detection  volumes.  In  the  small  detection  vol¬ 
ume  regime,  pi-OSCAR  results  in  a  —  l)  — >> 
27%  increase  in  the  measured  spin  signal  versus 
iOSCAR. 


15 


VI.  IMPROVEMENT  OF 
PERFORMANCE  OF  ESR  MRFM 
APPARATUS 

Improvement  of  force  detection  sensitivity  of 
our  MRFM  apparatus  is  an  ongoing  effort  in  our 
drive  towards  detecting  magnetic  resonance  sig¬ 
nal  from  a  single  nuclear  spin.  We  have  already 
demonstrated  EPR  MRFM  detection  with  signal 
sensitivity  corresponding  to  ~  1-2  electron  spins. 
Detection  of  a  single  nuclear  spin  will  require  im¬ 
provement  in  force  detection  sensitivity  by  three 
orders  of  magnitude.  Routes  toward  this  goal  in¬ 
clude  improvements  in  the  optical  interferomet¬ 
ric  cantilever  displacement  detection  and  reduc¬ 
tion  of  the  cantilever  thermal  noise.  Our  effort  in 
these  directions  is  described  below. 


A.  Optical  displacement  detection 


The  displacement  of  a  micromechanical  can¬ 
tilever  detected  by  fiber-optic  interferometry 
measures  the  interference  between  two  laser 
beams  one  of  which  is  reflected  off  the  end  of  the 
optical  fiber  and  a  second  reflected  back  into  the 
optical  fiber  by  the  reflective  surface  of  the  can¬ 
tilever.  The  visibility  V  characterizing  the  rela¬ 
tive  intensities  of  the  two  beams  is  given  by: 


V  = 


RfRc 

Rf  +  Rc 


(24) 


where  Rf  is  the  fraction  of  power  reflected  at  the 
end  of  the  fiber  and  Rc  is  the  fraction  of  power 
reflected  at  the  cantilever.  The  visibility  is  a  pa¬ 
rameter  of  critical  importance  due  to  the  fact  that 
the  shot  nose  spectral  density  S'shot  of  displace¬ 
ment  detection,  a  quantity  limiting  detection  sen¬ 
sitivity,  depends  on  visibility  as 


S'shot  00 


1 

V' 


(25) 


Thus,  increasing  V  improves  sensitivity.  We  have 
concentrated  significant  effort  to  improving  can¬ 
tilever  visibility.  By  using  various  engineering 
measures  including  uses  of  similar  materials  in 
apparatus  construction  to  minimize  thermal  drift 


FIG.  11:  MRFM  force  noise  as  a  function  of  signal 
modulation  frequency  recorded  for  two  values  of  can¬ 
tilever  visibility.  It  can  be  seen  that  the  onset  of  the 
shot  noise  shifts  from  ~  4  Hz  to  ~  40  Hz  as  visibil¬ 
ity  increases.  The  overall  reduction  of  the  thermal 
noise  floor  at  low  frequencies  from  ~  20  aN/\/Hz  to 
~  10  aN/\/Hz  is  due  to  improved  quality  factor  Q  of 
the  MRFM  cantilever. 


and  by  improving  procedures  for  doing  cantilever- 
fiber  alignment  we  were  able  to  increase  cantilever 
visibility  from  routinely  achieved  4%  to  70%.  The 
result  of  this  improvement  is  presented  in  Fig.  11 
showing  MRFM  force  noise  as  a  function  of  signal 
modulation  frequency  recorded  for  two  values  of 
cantilever  visibility.  It  can  be  seen  that  the  on¬ 
set  of  the  shot  noise  shifts  from  ~  4  Hz  to  ^  40 
Hz  as  visibility  increases.  The  overall  reduction 
of  the  thermal  noise  floor  at  low  frequencies  from 
~  20  aN/v^Hz  to  ~  10  aN/\/Hz  is  due  to  im¬ 
proved  quality  factor  Q  of  the  MRFM  cantilever 
that  will  be  discussed  further.  The  shift  of  the  on¬ 
set  of  the  short  noise  allows  modulation  of  MRFM 
signal  at  higher  frequencies  which  avoids  the  low 
frequency  1  //  noise  which  grows  as  the  tip  ap¬ 
proaches  the  sample. 


B.  Reduction  of  cantilever  thermal  noise  by 
increasing  quality  factor  Q 

The  ultimate  sensitivity  limit  for  mechanical 
force  detection  is  imposed  by  the  thermal  noise  of 
the  mechanical  cantilever  used  for  force  detection. 
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Frequency  Offset  (Hz) 

FIG.  12:  Measured  EPR  signal  spectral  density  us¬ 
ing  MRFM  iOSCAR  protocol  with  the  modulation 
frequency  of  75  Hz.  The  signal  lifetime  rm  =  1.24  s 
is  measured  from  the  signal  linewidth.  The  are  under 
the  signal  curve  is  proportional  to  the  number  of  spins 
N  participating  in  the  interaction.  The  observed  sig¬ 
nal  corresponds  to  N  =  13  spins.  The  signal  is  much 
larger  than  the  experimental  noise  floor. 


The  resulting  thermal  force  noise  Fth  (in  units  of 
N / \/Hz)  depends  on  the  cantilever  quality  factor 
Q  as: 

Fth  ~  75'  <26) 

Thus  the  thermal  noise  floor  of  the  MRFM  can¬ 
tilever  can  be  reduced  by  increasing  Q.  While 
there  are  multiple  mechanisms  that  degrade  Q, 
one  of  particular  importance  when  the  cantilever 
approaches  the  sample  is  cantilever  surface  con¬ 
tamination  and  electric  charge  trapping.  To  re¬ 
duce  the  influence  of  these  effects  we  have: 

•  Improved  the  process  of  gluing  a  probe  micro 
magnet  on  the  cantilever 

•  Changes  to  the  procedure  for  shaping  the 
probe  micro  magnet  with  the  Focused  Ion 
Beam  (FIB) 

•  Improved  methods  of  handling  and  storing  the 
cantilever 

•  Improved  the  design  of  the  cryogenic  appara¬ 
tus  to  minimize  contamination 

These  measures  were  developed  through  interac¬ 
tions  with  our  collaborators  at  UIUC.  The  result¬ 
ing  improvement  of  the  cantilever  thermal  force 


noise  by  a  factor  of  ~  2  is  shown  in  Fig.  11. 


C.  Verification  of  MRFM  apparatus 
performance 

Having  implemented  the  improvements  to 
MRFM  signal  detection  described  above,  we  have 
conducted  a  series  of  experiments  verifying  the 
performance  of  the  MRFM  apparatus.  These  ex¬ 
periments  were  also  necessitated  by  the  tempo¬ 
rary  shut-down  of  the  research  program  in  the 
first  half  of  2009  and  the  loss  of  the  personnel 
associated  with  the  project.  The  experiments 
were  conducted  on  the  standard  7-irradiated  sil¬ 
ica  sample  with  the  resulting  density  of  E'  cen¬ 
ters  of  6  x  1017  cm-3.  These  results  are  shown 
in  Fig.  12  which  presents  measured  EPR  signal 
spectral  density  using  MRFM  iOSCAR  protocol 
with  the  modulation  frequency  of  75  Hz.  The 
signal  lifetime  rm  =  1.24  s  is  measured  from  the 
signal  linewidth.  The  are  under  the  signal  curve 
is  proportional  to  the  number  of  spins  N  partic¬ 
ipating  in  the  interaction.  The  observed  signal 
corresponds  to  N  =  13  spins.  The  signal  is  much 
larger  than  the  experimental  noise  floor.  This  re¬ 
sult  demonstrates  force  detection  sensitivity  suffi¬ 
cient  for  our  research  program  including  END  OR 
MRFM  detection. 


VII.  PREPARATION  OF  CUSTOM 
N-IMPLANTED  DIAMOND  SAMPLES 
FOR  HIGH  SENSITIVITY/RESOLUTION 
ESR  MRFM  IMAGING 

At  present  there  are  only  two  large  manufac¬ 
tures  of  synthetic  single  crystal  bulk  diamond 
suitable  for  scientific  research:  Element  6  [25] 
and  Sumitomo  [26].  Both  companies  grow  dia¬ 
monds  via  either  a  high  pressure  and  high  tem¬ 
perature  (HPHT)  method  simulating  the  growth 
conditions  of  natural  diamond  but  at  a  faster  rate 
or  via  chemical  vapor  deposition  (CVD).  Chem¬ 
ical  vapor  deposition  method  of  growth  creates 
diamond  crystals  that  are  of  very  high  purity. 
Typically,  CVD  diamond  samples  from  these  two 
companies  contain  less  than  1  part  per  million  (1 
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FIG.  13:  The  results  of  the  theoretical  design  of  a 
custom  ion  implanted  diamond  sample  conducted  us¬ 
ing  publicly  available  SRIM  and  TRIM  software,  (a) 
Depth  of  nitrogen  ion  implantation  in  diamond  vs 
nitrogen  ion  energy.  Error  bars  represent  the  ex¬ 
pected  spread  of  the  implantation  depth,  (b)  The¬ 
oretically  designed  diamond  ion  implantation  density 
of  ^3  ppm  uniform  over  ~400  nm  ion  penetration 
depth  (shown  with  black  line)  achieved  by  combining 
shown  nitrogen  ion  implantation  profiles  at  multiple 
implantation  energies.  This  ion  implantation  scheme 
was  used  to  fabricate  the  custom  diamond  sample 
used  in  the  research  presented  in  this  report. 


ppm  =  2  x  10“ 17  cm~ 3)  of  nitrogen  PI  centers 
and  about  10-5  ppm  of  N-V  centers.  The  high 
pressure  and  high  temperature  diamond  contain 
significantly  more  nitrogen  impurities.  Using  the 
continuous  wave  EPR  we  measured  the  PI  center 
densities  of  the  HPHT  diamond  from  Element  6 
to  be  about  200  ppm  and  from  Sumitomo  to  be 
about  90  ppm.  Neither  company  guarantees  the 
density  of  PI  centers  in  delivered  HPHT  samples. 
As  a  result,  the  density  of  PI  centers  can  vary  sig¬ 
nificantly  in  a  commercial  HPHT  diamond  sam¬ 


ple  thus  making  such  a  sample  not  very  useful 
for  the  purposes  of  ESR  FMRFM  development. 
In  addition,  high  density  of  PI  centers  results  in 
fast  electron  spin  relaxation  thus  making  such  a 
system  a  poor  choice  as  a  sample  for  the  study  of 
electron  spin  dynamics. 

As  a  result,  diamond  samples  with  the  PI  cen¬ 
ter  density  in  the  range  of  1-100  ppm  cannot  be 
purchased  commercially.  However,  diamond  sam¬ 
ples  with  PI  center  density  of  ~10  ppm  would  be 
the  optimal  choice  for  the  study  of  electron  spin 
dynamics  since  they  possess  still  relatively  long 
spin  relaxation  time  but  can  already  exhibit  ef¬ 
fects  of  spin-spin  interaction.  In  response  to  this 
challenge,  we  have  decided  to  fabricate  custom 
diamond  samples  with  the  well  defined  PI  center 
density  by  custom  nitrogen  ion  implantation  of 
commercial  CVD  diamond  samples. 

Ion  implantation  at  a  single  ion  energy  results 
in  a  nonuniform  spatial  distribution  of  the  de¬ 
posited  atoms  (nitrogen  in  our  case),  while  we  are 
interested  in  a  uniform  spatial  distribution  of  ni¬ 
trogen  atoms  in  a  diamond  sample.  In  this  light 
we  have  theoretical  designed  a  custom  ion  im¬ 
planted  diamond  sample  using  publicly  available 
SRIM  and  TRIM  software  [27].  Fig.  13  shows 
the  resulting  theoretically  designed  diamond  ion 
implantation  density  profile  of  ^3  ppm  uniform 
over  ~400  nm  ion  penetration  depth  (shown  with 
black  line)  achieved  by  combining  shown  nitrogen 
ion  implantation  profiles  at  multiple  implantation 
energies.  This  ion  implantation  scheme  was  used 
to  fabricate  the  custom  diamond  sample  used  in 
the  research  presented  in  this  report. 

Nitrogen  ion  implantation  was  done  at  the 
commercial  implantation  service  Leonard  Kroko 
Implantation  [28]  following  the  implantation 
scheme  devised  using  SRIM  and  TRIM  software. 
During  the  implantation  the  sample  was  at  room 
temperature.  After  implantation  the  sample  was 
annealed  at  900  C  in  the  forming  gas  atmosphere 
(4%  H.2  in  Ar)  in  order  to  repair  damage  intro¬ 
duced  by  ion  irradiation  [29,  30].  Using  this  ap¬ 
proach  we  have  fabricated  two  types  of  test  sam¬ 
ples  for  ESR  MRFM 

•  Sample  with  uniform  ion  implantation  profile 

over  the  entire  sample  surface 

•  Samples  with  patterned  implantation  regions 
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fabricated  via  ion  implantation  through  a  gold 
mask  deposited  on  sample  surface 
Uniform  samples  allow  rapid  verification  of  the 
results  of  ion  implantation  since  they  do  not  re¬ 
quire  to  search  for  the  ion  implantation  region 
once  the  sample  is  cooled  to  liquid  He  temper¬ 
ature.  On  the  other  hand  patterned  samples 
present  a  much  more  interesting  target  for  ESR 
MRFM  imaging  and  for  the  study  of  dynamic 
spin  effects  such  as  spin  diffusion. 

The  results  of  the  experimental  studies  of  the 
custom  fabricated  samples  will  be  presented  in  a 
manuscript  ’’The  effect  of  spin  transport  on  life¬ 
time  in  nanoscale  systems“  that  is  currently  under 
review  in  Nature  Nanotechnology. 


VIII.  SIMULTANEOUS  MEASUREMENT 
OF  PI  AND  NV  CENTERS  IN  DIAMOND 

This  experiment  was  using  a  custom  fabricated 
uniformly  implanted  diamond  sample  with  the 
designed  PI  center  density  of  3  ppm.  It  is  known 
that  nitrogen  ion  implantation  followed  by  an¬ 
nealing  results  also  in  creation  of  N-V  centers 
[31].  Our  experimental  data  shown  in  Fig.  14 
demonstrates  simultaneous  detection  of  PI  and 
N-V  centers  in  diamond.  The  experimental  data 
(shown  with  green  dots)  clearly  exhibits  a  dou¬ 
ble  Lorentzian  structure  signifying  existence  of 
two  electron  spin  populations  with  different  life¬ 
times.  The  blue  and  red  curves  are  the  individual 
Lorentzians  fitting  the  experimental  curve.  The 
blue  curve  has  a  longer  lifetime  of  220  ms  and 
corresponds  to  a  fewer  spins  which  is  the  signa¬ 
ture  of  N-V  center  population.  The  red  curve, 
corresponding  to  a  larger  spin  population,  has  a 
shorter  lifetime  of  12  ms  and,  undoubtedly,  orig¬ 
inates  from  PI  centers.  Based  on  this  data,  we 
estimate  that  the  density  of  the  observed  N-V 
centers  is  ^0.1  ppm  which  is  in  a  good  agree¬ 
ment  with  the  expected  N-V  center  generation 
rate  [31].  This  experimental  data  demonstrates 
the  ability  of  ESR  MRFM  to  simultaneously  de¬ 
tect  and  distinguish  between  two  electron  spin 
populations  with  different  lifetimes.  Such  a  capa¬ 
bility  is  important  for  the  study  of  spin  systems 
of  practical  scientific  significance. 


FIG.  14:  A  force  power  spectrum  density  (PSD)  on 
the  MRFM  cantilever  demonstrating  simultaneous 
detection  of  PI  and  N-V  centers  in  diamond.  The 
experimental  data  (shown  with  green  dots)  clearly 
exhibits  a  double  Lorentzian  structure  signifying  ex¬ 
istence  of  two  electron  spin  populations  with  different 
lifetimes.  The  blue  and  red  curves  are  the  individual 
Lorentzians  fitting  the  experimental  curve.  The  blue 
curve  has  a  longer  lifetime  of  220  ms  and  corresponds 
to  a  fewer  spins  which  is  the  signature  of  N-V  center 
population.  The  red  curve  corresponding  to  a  larger 
spin  population  has  a  shorter  lifetime  of  12  ms  and, 
undoubtedly,  originates  from  PI  centers.  The  blue 
curve  is  offset  for  clarity. 


IX.  ENDOR  AND  EPR 
CHARACTERIZATION  OF  Pi  CENTERS 
IN  DIAMOND 

Electron  Nuclear  Double  Resonance  (ENDOR) 
is  a  very  attractive  approach  to  measuring  an  in¬ 
dividual  nuclear  spin  magnetic  resonance  through 
detecting  the  dynamics  of  an  electron  spin  cou¬ 
pled  to  the  nuclear  spin  through  hyperfine  cou- 
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EPR  and  NMR  Transitions 


FIG.  15:  Energy  level  diagram  of  the  14 N  Pi  center  in 
diamond  with  the  electron  spin  of  1/2  and  the  nuclear 
spin  of  1. 


pling.  In  this  method,  the  requirements  on  the 
MRFM  force  detection  sensitivity  is  relaxed  due 
to  the  fact  the  electron  magnetic  moment  fie  is 
approximately  1000  times  greater  than  that  of 
an  individual  nuclear  spin  jin.  Though  this  type 
of  nuclear  spin  detection  will  be  applicable  only 
to  systems  with  strong  electron-nuclear  hyper- 
fine  coupling,  demonstration  of  an  individual  nu¬ 
clear  spin  detection  using  ENDOR  will  be  an  im¬ 
portant  step  towards  detecting  individual  nuclear 
spin  directly. 

We  have  chosen  nitrogen  doped  diamond  as  a 
test  system  for  demonstrating  ENDOR  nuclear 
spin  detection  using  MRFM.  Diamond  is  com¬ 
posed  of  two  interpenetrating  face  centered  cu¬ 
bic  (FCC)  structures  of  carbon  atoms.  Doping 
a  diamond  crystal  with  nitrogen  results  in  the 
formation  of  so-called  Pi  centers  in  which  a  ni¬ 
trogen  atom  substitutes  for  carbon.  Substitution 
of  Nitrogen  with  its  extra  valence  electron  rel¬ 
ative  to  carbon  results  in  an  unpaired  electron 
shared  by  the  nitrogen  and  one  of  the  four  carbon 
atoms.  Hence  there  are  four  equivalent  Pi  centers 
aligned  along  one  of  four  directions:  (111),  (111), 
(111),  and  (111).  The  resulting  energy  configu¬ 
ration  of  the  Pi  center  will  depend  on  whether 
the  nitrogen  atom  is  147V  with  nuclear  spin  1=1 
or  15N  with  nuclear  spin  I  =  1/2.  Since  14N 
has  natural  abundance  of  99.632%  the  most  com¬ 
mon  electron-nuclear  configuration  of  the  Pi  cen¬ 
ter  has  5  =  1/2  and  1=1  with  the  resulting 
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FIG.  16:  Conventional  CW  Pi  center  EPR  in  dia¬ 
mond  measured  as  a  function  of  the  sample  orien¬ 
tation  relative  the  applied  magnetic  field  measured 
at  temperature  T  =  300  K  and  microwave  frequency 
/  =  9.4  GHz.  The  orientation  of  the  applied  magnetic 
field  was  varied  relative  to  the  (100)  axis  of  the  dia¬ 
mond  sample.  When  magnetic  field  is  aligned  along 
(100),  all  four  types  of  Pi  centers  form  the  same  angle 
P  with  the  magnetic  field  thus  resulting  three  identi¬ 
cal  EPR  lines.  As  magnetic  field  is  rotated  towards 
the  (010)  direction,  Pi  centers  split  in  two  groups  by 
the  magnitude  of  the  resulting  angle  P  which  is  man¬ 
ifested  by  splitting  Iz  =  ±1  EPR  transitions  into  two 
lines. 


electron-nuclear  Hamiltonian  H  given  by: 

=  veSz  +  SZIZ^A\  cos2  /3  +  sin2  ~/3 
+PI2Z  cos2  /3  -vnIz, 

where  ve  is  the  electron  spin  resonance  frequency, 
A\\  and  A±_  are  the  hyperfine  coupling  parallel 
and  perpendicular  to  the  direction  of  the  carbon- 
nitrogen  bond,  P  is  the  nuclear  quadrupole  cou¬ 
pling,  vn  is  the  nitrogen  nuclear  magnetic  reso¬ 
nance  frequency,  P  is  the  angle  made  by  the  ap¬ 
plied  magnetic  field  with  the  direction  of  carbon- 
nitrogen  bond,  Sz  and  Iz  are  the  z  (parallel  to 
applied  magnetic  field)  components  of  the  elec¬ 
tron  and  nuclear  spin.  The  resulting  Pi  energy 
diagram  is  shown  in  Fig.  15.  These  energy  char¬ 
acteristics  make  Pi  centers  in  diamond  an  ideal 
sample  for  demonstrating  ENDOR  detection  with 
MRFM. 

We  have  acquired  a  nitrogen  doped  diamond 
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FIG.  17:  Observed  splitting  of  Iz  —  ±1  peaks  as  a 
function  of  orientation  of  applied  magnetic  field.  Blue 
squares  represent  measured  peak  splitting  and  solid 
line  is  the  fit  of  the  data  with  Eq.  27  using  A\\  and 
A±  as  fitting  parameters.  The  independent  variable 
is  the  angle  made  by  applied  field  Ho  with  (100)  axis. 


sample  from  a  commercial  supplier,  Element  Six 
(http://www.e6.com).  The  sample  has  dimen¬ 
sions  2.7  mm  x  2.7  mm  x  0.5  mm  with  an  ex¬ 
pected  Pi  center  concentration  of  200  ppm.  Prior 
to  the  MRFM  experiment  it  is  essential  to  charac¬ 
terize  the  sample  by  conventional  EPR/ENDOR 
means  with  the  main  objective  to  verify  a  suffi¬ 
cient  concentration  of  Pi  centers  and  confirm  the 
strength  of  the  hyperfine  coupling.  These  experi¬ 
ments  were  conducted  using  Bruker  Electron  Spin 
Resonance  (ESR)  Spectrometer  acquired  with 
funds  from  the  ARO  DURIP  grant  W911NF-07- 
1-0305.  The  experimental  results  are  presented 
below. 

With  availability  of  four  equivalent  Pi  cen¬ 
ters  indicated  by  three  allowed  EPR  transitions 
each  (See  Fig.  15),  the  resulting  EPR  spectrum  is 
strongly  dependent  on  the  orientation  of  the  ap¬ 
plied  external  magnetic  field  relative  to  the  axis 
of  the  diamond  crystal.  Depending  on  the  ori¬ 
entation  of  the  field  the  EPR  transitions  corre¬ 
sponding  to  Iz  =  ±1  can  be  split  by  the  hyper¬ 
fine  coupling  into  up  to  four  EPR  lines  due  to 
four  different  angles  between  the  magnetic  field 
and  the  direction  of  the  carbon-nitrogen  bond. 
This  dependence  can  be  exploited  to  measure  the 
strength  of  the  hyperfine  coupling  as  shown  in 
Fig.  16  showing  CW  EPR  measurements  done 
on  Pi  centers  diamond  at  temperature  T  =  300 
K  and  microwave  frequency  /  =  9.4  GHz.  In 
the  experiment  EPR  spectra  were  measured  as  a 
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FIG.  18:  Conventional  EPR  and  ENDOR  spectra  of 
Pi  center  in  diamond  acquired  with  applied  magnetic 
field  Ho  aligned  along  the  (100)  direction  of  the  dia¬ 
mond  crystal,  a)  Conventional  Pi  EPR  signal  demon¬ 
strating  EPR  transitions  indicated  in  Fig.  15  acquired 
at  applied  /x-wave  frequency  f^w  =  9.4  GHz.  b) 
Conventional  Pi  ENDOR  spectra  acquired  with  ex¬ 
ternal  magnetic  field  ifo,  combined  with  the  /x-wave 
frequency  f^w  =  9.4  GHz,  set  to  select  one  of  the 
EPR  transitions  with  Iz  =1,  0  or  -1  (as  indicated  in 
Fig.  15).  ENDOR  spectra  are  shown  as  a  function 
of  frequency  of  applied  rf  radiation.  ENDOR  spec¬ 
tra  a  offset  for  clarity.  The  resulting  spectra  have 
all  the  features  expected  from  the  hyperfine  coupled 
electron-nuclear  system.  As  f^w  is  aligned  with  the 
side  peaks  (Iz  =  ±1)  of  the  EPR  spectrum,  we  ob¬ 
serve  signatures  of  two  NMR  transitions  (as  can  be 
seen  in  Fig.  15).  Each  transition  is  split  into  four  in¬ 
dividual  lines  due  to  slight  misalignment  of  magnetic 
field  with  the  (100)  direction  resulting  four  slightly 
different  angles  with  carbon-nitrogen  bond.  As  f^w 
is  aligned  with  the  central  (Iz  =  0)  EPR  peak,  we 
observe  four  NMR  transitions  (as  can  be  seen  in 
Fig.  15).  Each  transition  is  yet  again  split  into  four 
lines  due  to  field  misalignment.  In  the  ENDOR  spec¬ 
trum  shown  in  b)  some  NMR  lines  overlap  thus  some¬ 
what  obscuring  the  observed  picture. 


function  of  the  sample  orientation  relative  the  ap¬ 
plied  magnetic  field.  As  magnetic  field  is  aligned 
with  the  (100)  direction  of  the  crystal,  all  four 
types  of  Pi  centers  form  the  same  angle  (3  with 
the  magnetic  field  thus  resulting  three  identical 
EPR  lines.  As  magnetic  field  is  rotated  towards 
the  (010)  direction,  Pi  centers  split  in  two  groups 
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by  the  magnitude  of  the  resulting  angle  /3  which  is 
manifested  by  splitting  Iz  =  ±1  EPR  transitions 
into  two  lines  as  demonstrated  in  Fig.  16. 

The  angular  dependence  of  the  Iz  =  - bl  line 
splitting  A K  is  described  by: 

A K  =  yj A ^  cos2  pi  +  A\  sin2  /?i 

-  yj |  cos2  +  A2±  sin2  /?2 ,  (27) 

where  and  /?2  are  the  angles  between  the  mag¬ 
netic  field  and  the  carbon-nitrogen  bond  for  the 
two  groups  of  P\  centers.  This  equation  was  used 
to  fit  experiments  data  as  shown  in  Fig.  17  show¬ 
ing  observed  splitting  of  Iz  db  1  peaks  as  a  func¬ 
tion  of  orientation  of  applied  magnetic  field.  Blue 
squares  represent  measured  peak  splitting  and 
solid  line  is  the  fit  of  the  data  with  Eq.  27  us¬ 
ing  A\\  and  A±  as  fitting  parameters.  The  angle 
on  the  x-axis  is  the  angle  made  by  applied  field 
H0  with  (100)  axis. 

The  values  of  longitudinal  and  transverse  hy- 
perfine  splitting  obtained  from  this  data  fit  are 
A\\  =  40.7  ±  0.1  G  and  A±  =  29.3  ±  0.1  G. 
These  results  are  in  a  good  agreement  with  the 
values  available  in  literature  [32]  with  A\\  =  40.24 
G,  A±=  29.73  G.  These  results  also  confirm  that 
the  sample  has  a  sufficient  Pi  center  density  and 
is  therefore  suitable  for  END  OR  MRFM  experi¬ 
ments. 

We  have  conducted  conventional  END  OR  ex¬ 
periments  using  our  CW  EPR  spectrometer;  re¬ 
sults  are  shown  in  Fig.  18.  The  applied  magnetic 
field  H0  was  aligned  parallel  to  (100)  resulting 
in  a  simple  EPR  spectrum  with  three  EPR  lines 
(Fig.  18a).  The  microwave  frequency  f^w  was  set 
to  a  fixed  value  matching  one  of  the  three  EPR 
transitions.  With  f^w  fixed,  an  ENDOR  spec¬ 
trum  was  recorded  as  a  function  of  the  frequency 
fr f  of  the  applied  rf  radiation  generated  by  an 
additional  rf  source.  The  resulting  spectra  have 
all  the  features  expected  from  the  hyperfine  cou¬ 
pled  electron-nuclear  system.  As  f^w  is  aligned 
with  the  side  peaks  (Iz  =  ±1)  of  the  EPR  spec¬ 
trum,  we  observe  signatures  of  two  NMR  transi¬ 
tions  (see  Fig.  15).  Each  transition  is  split  into 
four  individual  lines  due  to  a  slight  misalignment 
of  the  magnetic  field  with  the  (100)  direction  re- 


FIG.  19:  LAO/STO  sample  hysteresis  curve,  shown 
with  diamagnetic  background  subtracted,  measured 
using  a  SQUID  magnetometer.  The  data  was  ac¬ 
quired  at  300  K. 


suiting  four  slightly  different  angles  with  carbon- 
nitrogen  bond.  When  f^w  matches  the  central 
(Iz  =  0)  EPR  peak,  we  observe  four  NMR  transi¬ 
tions  (as  can  be  seen  in  Fig.  15).  Each  transition 
is  yet  again  split  into  four  lines  due  to  field  mis¬ 
alignment,  hence  we  see  16  transitions  in  all.  In 
the  ENDOR  spectrum  shown  in  Fig.  18b  some 
NMR  lines  overlap  thus  somewhat  obscuring  the 
picture. 

In  conclusion,  we  have  successfully  character¬ 
ized  nitrogen  doped  diamond  sample  and  have 
confirmed  its  suitability  for  use  in  ENDOR 
MRFM  experiments. 


X.  INTERFACE  MAGNETISM  IN 
LAO/STO 

We  have  conducted  preliminary  investigations 
of  physical  systems  that  might  serve  as  sam¬ 
ples  for  future  ESR-MRFM  studies.  The  sys¬ 
tems  of  particular  interest  would  exhibit  mag¬ 
netic  properties  at  nanoscale.  Such  a  system 
would  require  use  of  ESR-MRFM  for  its  stud¬ 
ies,  since  the  conventional  magnetic  resonance 
techniques  do  not  have  sensitivity  necessary  for 
the  study  of  nanoscale  samples,  thus  highlight¬ 
ing  the  advantages  of  the  ESR-MRFM  over  other 
experimental  approaches.  One  of  the  promis- 
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ing  systems  in  this  regard  is  the  2D  interface 
in  LaA10s/SrTi03(LA0/ST0)  system  where  in¬ 
teresting  magnetic  properties  have  been  reported 
including  coexistence  of  ferromagnetism  and  su¬ 
perconductivity  [33].  High  spin  sensitivity  and 
the  high  spatial  resolution  of  MRFM  will  be  in¬ 
valuable  for  the  study  of  such  a  system.  As  it 
is  common  with  with  ESR-MRFM  applications, 
the  potential  samples  have  to  be  characterized  by 
conventional  techniques. 

We  performed  a  series  of  SQUID  measurements 
on  an  LAO/STO  sample.  The  sample  was  grown 
by  the  group  of  Prof.  C.  H.  Ahn  at  Yale  Univer¬ 
sity  (New  Haven,  CT).  The  sample  has  metallic 
interface  states  that  were  confirmed  by  measur¬ 
ing  the  conductivity  as  a  function  of  temperature, 
finding  a  sheet  resistance  of  1000  D/square  at  a 
temperature  of  2  K. 

We  collected  hysteresis  curves  to  investigate 
the  ferromagnetism  as  reported  in  the  literature 
[34].  For  our  measurements,  we  took  both  field 


cooled,  and  zero-field  cooled  scans  at  a  variety 
of  temperatures,  since  the  literature  reports  an 
increase  in  magnetic  moment  at  low  tempera¬ 
tures  with  the  example  of  acquired  hysteresis 
curve  shown  in  Fig.  19.  Our  measurements  show 
a  strong  diamagnetic  background,  which  we  at¬ 
tribute  to  the  bulk  STO  substrate,  but  they  also 
reveal  weak  ferromagnetism  that  is  roughly  inde¬ 
pendent  of  temperature.  This  result  is  surpris¬ 
ing  given  that  the  moment  is  expected  to  grow 
with  decreasing  temperature.  Upon  further  in¬ 
vestigation,  we  concluded  that  this  effect  must 
be  due  to  contaminants  on  the  unpolished  sur¬ 
faces  of  the  STO  substrate,  which  has  been  found 
to  yield  weak  ferromagnetic  signal  [35].  We  can 
conclude  that  that  the  net  magnetic  moment  at 
the  LAO/STO  interface  was  smaller  than  10-6 
emu,  the  limit  of  the  SQUID  that  was  used  for 
measurement.  Thus,  this  might  mean  that  ESR- 
MRFM  might  be  one  of  the  few  methods  for  the 
study  of  magnetism  at  LAO/STO  interface. 
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